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ABSTRACT 

A systematic study of the synchrotron X-ray emission from supernova remnants (SNRs) has been conducted. 
We selected a total of 12 SNRs whose synchrotron X-ray spectral parameters are available in the literature 
with reasonable accuracy, and studied how their luminosities change as a function of radius. It is found that 
the synchrotron X-ray luminosity tends to drop especially when the SNRs become larger than ^^5 pc, despite 
large scatter This may be explained by the change of spectral shape caused by the decrease of the synchrotron 
roll-off energy. A simple evolutionary model of the X-ray luminosity is proposed and is found to reproduce the 
observed data approximately, with reasonable model parameters. According to the model, the total energy of 
accelerated electrons is estimated to be lO"^^""*** ergs, which is well below the supernova explosion energy. The 
maximum energies of accelerated electrons and protons are also discussed. 
Subject headings: acceleration of particles — ISM: supernova remnants — X-rays: ISM 
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1. INTRODUCTION 

Young supernova remnants (SNRs) are widely beli eved to 
be the main source of Galactic cosmic rays. Kovam Tet al.l 
(11995) discovered synchrotron X-rays from shells of 
SN 1006, which was the first observational clue of cosmic- 
ray electrons being accelerated up to the TeV energy range. 
Later, several young SNRs were found t o have synchrotron 
X-ray shells (c.f QCoyama et al. 1997t iSlane etall IT999t 
iBamba et al]|2000t mnk & Laming 2003"). Another piece of 
evidence for cosmic-ray acceleration in SNRs was obtained 
from gamma-ray observations. Very-high-energy (VHE) and 
GeV gamma-rays have been detected from several SNRs, al- 
though, it is still unclear whether their origin is hadronic 
or leptonic_(e.g., lAharonian et all 12004 l2007t lAbdo et all 
UOlOb c, 26m. 

Despite these various pieces of evidence in favor of acceler- 
ation, it is still unclear how the acceleration process evolves in 
SNRs, in particular how the accelerated particles cool down 
and how they escape from SNRs. One reason for this lack of 
understanding is that past observational studies have concen- 
trated on individual sources alone. In this paper, we investi- 
gate, for the first time, the time evolution of SNR synchrotron 
X-rays, using all available data in literature for sources in 
which this emission component is bright. In section |2] we 
describe the evolution of synchrotron X-ray luminosity in our 
sample. An interpretation and a simple model to support it are 
in section[3] Finally, section|4]is devoted to discussion of the 
results. 

2. SAMPLE SELECTION AND RESULTS 

We searched the literature for reports of synchrotron X- 
ray emission from SNRs and found a total of 14 such 
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sources. However, we could not use all of these for the cur- 
rent study, for variou s reasons. Kepler ("Bamba et al.|j2005t 
[Reynolds etalll2007h and G330.2H-1.0 (Park et al. 200l) had 
to be removed from our sample, because of the difficulty of 
estimating the total synchrotron X-ray luminosity, caused by 
a low detection significance or by contamination from thermal 
X-rays. Thus our sample consists of the 12 SNRs listed in Ta- 
ble [J We used the latest value of the 2-10 keV synchrotron 
X-ray unabsorbed flux from the references listed in the table. 
When the reported flux was not absorption-corrected or was 
for a different energy band, we normalized it unabsorbed flux 
in the 2-10 keV using the best-fit model in the references. 
Synchrotron X-ray emission may have a spectral roll-off in 
the 2-10 keV band. When significant roll-off was reported in 
the references, we took such a component into account while 
computing the flux. Some SNRs (Cas A, Tycho, SN 1006) 
have bright thermal X-ray emission, which makes it difficult 
to estimate the synchrotron X-ray flux. In these cases, we 
inferred the synchrotron X-ray flux based on wide-band in- 
formation and/or detailed spectroscopy to isolate synchrotron 
X-rays from thermal component, as we believe such analysis 
can provide the most reliable results currently available. Two 
of the sample sources, W28 and G156.7H-5.7, are largely ex- 
tended but were only partly observed with enough exposure 
in X-rays. As the X-ray luminosities can be deduced only for 
the observed regions (5% for W28 and 26% for G156.7H-5.7, 
respectively), possible contributions from the unobserved re- 
gions are included in the upper error bars. When we estimated 
these contributions, we assumed that the surface brightness of 
the unobserved region in each source is identical to that of 
the observed one. The distance uncertainty is also included 
into the errors on the luminosity. We used a 10% distance un- 
certainty for cases where no uncertainty was available in the 
literature. In summary, we took into account the following 
uncertainties; statistical errors given in literatures, systematic 
distance uncertainty, and the limited coverage of the SNR ex- 
tension. 

In order to study the time evolution, we need to know the 
SNR ages. However, we have only four SNRs whose ages are 
historically known. Another parameter, the ionization time 
scale of heated plasma, is also difficult to determine since 
several SNRs show no thermal emission. We thus use the 
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physical radi us of e a ch SN R, R^, as the age indicator 
is taken from iGreenI ( 120091) . The distance uncertainty is in- 
cluded in the error computation. Some SNRs have distorted 
shapes, as shown in Tab. [T] This is also included in the errors 
on the radii. All the derived parameters are shown in Tab. [1] 
The radius of a SNR depends on the density of the interstellar 
medium, but rather insensitively. For remnants in the Sedov 
phase, Rs oc no ''^^ (where no is the upstream number density). 
The radius changes by only a factor of ^2 even if the up- 
stream density changes by 2 orders of magnitude. The radius 
also has very weak dependency on the explosion energy Eq, 

Rs oc £'0''^. On the other hand, Rs oc fage^''^ (where f^ge is the 
age of the remnant). Thus the radius is a good age indicator 

Figure[T]shows the synchrotron X-ray luminosity as a func- 
tion of the radius. One can see that the luminosity is in 
the order of lO^'^-lO^^ ergs s"' when the SNRs are smaller 
than Rs ^5 pc, which corresponds to an age of a few hun- 
dred years, whereas it decreses to 10^^-10^^ ergs s"' beyond 
Rs ^5 pc. it appears that the luminosity drops off rapidly at 
Rs ^5 pc, although the scatter is rather large. Note that all 
the SNRs with Rs < 5 pc display synchrotron X-rays brighter 
than lO^'* ergs s~' (lGreenl20 09). whereas most of larger SNRs 
do not have significant synchrotron X-rays. This fact makes 
the difference of luminosities larger between the two regions. 
The drop off in non-thermal X-ray luminosity reaches two 
or three orders of magnitude, which is much larger than 
the errors on the individual luminosities. We have carried 
out a similar analysis for the thermal X-ray luminosity for 
SNRs in the Large Magellanic Cloud (Williams et al. 1999) 
and the rad io luminosity of Galactic SNRs (Green i2009l: 
iCase & Bhattacharva 1998.). However, no drop off like that 
observed for synchrotron X-rays has been found. 

3. EVOLUTION OF SYNCHROTRON X-RAYS 

Here we consider the cause of the decrease in synchrotron 
X-ray emission identified in the previous section. The i^/V" 
spectrum of synchrotron X-ray radiation has a peak around 
the roll-off frequency, i^ion, above which the flux rapidly de- 
cays towards higher fre quencies. As SNRs evolve, Vi-ou is 
known to decrease (e.g.. iBamba et al.l 12005). The observed 
rapid decay in synchrotron X-rays around Rs ^ 10 pc could 
be caused by i/roii passing through the X-ray band to lower en- 
ergies. In this section, we discuss this possibility in detail. We 
introduce essential argument first, and show a simple model 
to support it later 

We assume that the electron acceleration is energy-loss- 
limited, in which the maximum energy of electrons, /Sg™", 
is determined from the balance of the synchrotron loss and 
acceleration: 

^■"'"=1^(10^) (life)"''' <" 

where ^ is the gyro-factor In deriving Eq. ([T]l, we 
equate the acceleration time, facc(£) = 2Q^cE /3eVs^Bd 
with the synchrotron cooling time, fsyn(£) = 
125 yr(£/10 TeV)-i(Brf/100 fiGT^, where v, and 
are the shock velocity and the downstream magnetic field, 
respectively. In this case, we derive 

/it'roii ~0.4keVr'(v,/10^cms-')2 , (2) 

which is independent of Bd (e.g..' Aharonian & Atovanll 19991: 
lYamazaki et al...2006.) . The roll-off frequency and the shock 



velocity ha ve b een measured in several SNR s (Cas A, 
SN 1 006; [Reyn olds & Keohane 1999; Patnau de & FesenI 
120091; iBambaet al.»2008i; iGhavamian et al.,.2002,) . For these 
SNRs, eq.dljl is consistent with the observational values to 
within 1 order of magnitude assuming ^ = 1 . Also note that the 
effect of par ticle escape from the shock region (e.g.. Reyn olds! 
1998; Ohir a et al.ll2010l) is not considered in this paper. Thi s 
effect might be important for older SNRs dOhira et al .11201 ll) . 
However, our present model reproduces the observed trend 
well, suggesting that particle escape is not yet significant for 
these SNRs (see also the last paragraph in section 4). 

The X-ray spectrum of synchrotron radiation is well ap- 
proximated analytically. It is mainly determined by shock 
dynamics in the SNR if the synchrotron X-ray-emitting elec- 
trons with energy E satisfy fsyn(£^) < fage, where fage is the age 
of the SNR. In this case, the energy spectrum of electrons is 
generally given by 

N,{E)=AE-P(l+E/Eh)-'exp[-{E/Er'"f] , (3) 

where the break energy, £/,, is determined by fsyn(£&) = fage. 
During the acceleration, electrons with E > Et suffer en- 
ergy loss via synchrotron cooling, which causes a steepen- 
ing of the energy spectrum (e.g., ILongaii1ll994i) . We also 
note that the shape of the cutoff is analytically obtained in 
IZirakashvili & AharonianI (l2007h . Given the electron distri- 
bution, we calculate the approximate formula of the X-ray 
luminosity, L^. The characteristic frequency, v^y^iE), of the 
synchrotron radiation emitted by electrons with energy E is 
given by 

/!z^syn(£)--0.12keV(Bj/10/iG)(£/10TeV)2 . (4) 

Then, as long ?&v<Vh = v^yn{Eb), we can apply the standard 
formula, o^AB^^^^^'^y-^P-^'^l'^ (e.g.. lLongaiij[T994,) . On the 
other hand, if Vb < i', the spectral slope steepens (L^ oc v'pI^) 
due to the steepening of the electron distribution, and we de- 
rive 

Uo,AB/P^'y^Vb-^P-'y\vlvbrPl^cxp(-./VJ^,) 
^AB/P-^'^l^v-Pl^ txp{-^v/v,,n) , (5) 

where we assume tha t Ly is continuous at v = Vh, a nd we again 
use the result of Zira kashvili & AharonianI (l2007h for the cut- 
off shape. Calculating vi, as 

0.19 ke V(Bd/ 10 AiG)-3 (fage/ lOVr^ , (6) 

we find that throughout the evolution of the SNR, the X- 
ray band (2-10 keV) always Ues above vi, because for young 
SNRs (fage ^10^ yr), the magnetic field may be amplified to 
Bd > 10 ^G re.g..'Bamba et al.' 2003. 2005VVin k & LamingI 
120031) . Equation. ^ is thus a good approximation of the X- 
ray luminosity for any arbitrary epoch. In particular, when 
p « 2, the X-ray luminosity is insensitive to the magnetic 
field, depending instead mainly on i/foii. As the SNR ages, the 
shock velocity Vj decreases and v^oW becomes smaller. Equa- 
tion (|2| tells us that v^ow is below the X-ray band (2-10 keV) 
when Vj < lO^cm s"', so that the X-ray luminosity drops off. 

In order to demonstrate the above argument, we construct 
a simple model to calculate the synchrotron X-ray flux. In 
our model, a simple shock dynamics scenario is considered. 
We assume that the forward shock velocity of SNRs Vs is a 
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function of the age of SNR tgge as follows: 

' V, (tage < h ; Free expansion phase) 

ft 

Vi[^] (fi < tage < Sedov phase) 



/, (7) 

Vil't] ( 1f ) te < fflgf ; Radiative phase), 



fi=2.1 xlo2(£5i/no)'/'v,-,9"'/' 
f2=4xl04£5i''''W/'' yr, 



yr, 



(8) 
(9) 



where v,- = v/g x 10 cm s" , ZSgj = £^5i x 10 erg, and no 
are the initial velocity, initial energy of ejecta, and the up- 
stream number d ensity, respectively (e.g. jBlondin et al.ll998l; 
[Vamazaki et al. 2006). The model assumes that the SNR for- 
ward shock propagates into a homogeneous medium. This 
is, however, not always in the case, since the SNRs could be 
located inside of superbubbles made by pre-explosion stellar 
winds. More detailed m odels for core-c ollapsed type SNRs, 
such as RCW 86 (e.g., Vink et al.|[T997b . should be appUed. 
But it is left for future work. 

As illustrated above, the X-ray luminosity has little depen- 
dence on the magnetic field. Nevertheless, we model the 
evolution of the downstream magnetic field, Bd, as follows. 
Usually, the amplified magnetic field Z?amp is simply assumed 
to scale with Vj (Volk et al. 2005; Vink 2008), because mag- 
netic field evolution remains ill-understood. Here we adopt 
an assumption similar to Vol k et al. (2005), leading to a de- 
pendence of Bj oc v,. Let the energy density of the amplified 
magnetic field Ub = Bamp^/87r be proportional to the thermal 
energy density U,h, 

■Samp . , 3 

— — = esUih = -eBrni)kTd , 
on Z 

where e^, r, and kTj are the energy-partition parameter, a 
compression ratio, and a downstream temperature obtained 
by the Rankine-Hugoniot equation kT^ = [(r- 1)/ r^]mpVs^, re- 
spectively. Then, we obtain Ba^p = [l2iTmp{r- l)eBMo/'"]'^^^.s- 
The downstream magnetic field is thus given by Bj = 
max{Batnp, rBisM}, where Bism = 10 /iG is the field strength 
in the interstellar matter. 

For given fage, we calculate according to Eq. (|7]), the 
shock radius /?, by the integration of v,, and the luminosity 
of synchrotron X-rays numerically by using the value of Bd 
and the electron distribution given by Eqs. ^ and (|3]l. We 
adopt £51 = V, 9 = 1, p = 2.0, r = 4, ^ = 1, and = 0.01 as fidu- 
cial parameters. Note that for ^ 0.01, the evolution of the 
amplified magnetic field in young SNRs is approximately re- 
produced (IVolk et al.ll2005h . This is also co nsistent with pre- 
vious observational implications for SNRs ( Ghavamian et al.l 
12002; Bam ba etal.l l2003; Parizotet al. 2006). One finds that 
with the above assumption for the magnetic field, the cool- 
ing time of synchrotron X-ray-emitting electrons is always 
smaller than the SNR age, faee, so that the X-ray luminosity 
is well approximated by Eq. 1©. Also note that the value of 
^ is near unity, implying that acceleration is near the Bohm 
limit. We assume that the normalization factor of the electron 
distribution A is constant with time. This is justified if the 
amount of accelerated particles is proportional to the prod- 
uct of the fluid ram pressure and the SNR volume, that is, 
A oc {n(,Vs^)Rj'. In the Sedov phase {Rs oc f^^^) during which 
frou crosses the X-ray band, we find that A is constant with 
time. 



Fig.[T]shows the results for no = 5.0, 1.0, and 0.1 cm"^. For 
the first two cases, we set A = 9.0 x 10'*^(no/l cm"^), whereas 
A = 9.0 X 10"^^ for no = 0.1 cm"^. A is expressed in the cgs 
units. One can see that our model roughly reproduces the 
observed trend with reasonable model parameters. The as- 
sumed density, no = 0.1-5.0 cm"-', is typical for the interstel- 
lar medium around SNRs, although there are discrepancies 
for some individual samples. Some SNRs are likely to be lo- 
cated in bu bbles or very low density regions. For example, 
IVink et al.i (f!997 ;) showed that RCW 86 is in a super-bubble 
and expands into extremely low density materials. As a re- 
sult, this SNR should have a very large radius compared with 
its age. A is also a free parameter which may change from 
source to source. Adjusting A to within 1 order of magni- 
tude, it is possible to find a solution for the observed density 
in each case. Since our aim is simply to reproduce the overall 
evolutionary trend in an approximate manner, we ignore such 
fluctuations in this work. 

4. DISCUSSION 

We have reproduced the observed fast drop-off in syn- 
chrotron X-ray lumin osity with SNR radius using a quite sim- 
ple model. Recently, IPatnaude et al.l ( 1201 lb discovered a de- 
cline of a few percent in the synchrotron X-ray power from 
one source, Cas A, and reached a similar conclusion. 

The total and maximum energy of accelerated particles can 
be derived from our model. With the adopted values of A 
and Eq. (|3]l, we can calculate the total energy of accelerated 
electrons with the energy above m^c^ to be 10'*'-10'*'* ergs, 
which is reasonable since this range is much smaller than 
Esy Note that these estimated values could increase by 1 
order of magnitude if we increase ^ by up to 10. Even if 
we alter the value of p, the results remain similar if we ap- 
propriatelyreset the normalization A. This is expected from 
equation ©. For example, the results for p = 2.2 are al- 
most perfectly in agreement with the case of p = 2.0 if we 
choose A = 1.7 X 10'*^(no/l cm~^) no = 5.0 and 1.0 cm~^ and 
A = 1.7 X 10"*^ for no = 0.1 cm"-'. In these cases, the total en- 
ergy of accelera ted electrons w ith energy above is again 
lO^^-io^s ergs. iBamba et all (|2003) estimated the energy of 
injected accelerated electrons in SN 1006 for the first time. 
This was found to be ^ 10'*'*-10'*^ ergs per small segments 
of filaments, although with rather large uncertainty. The to- 
tal energy of accelerated electrons is approximately 10"^^- 
10"*^ ergs considering the size of these segm ents, which is 
consistent with our result. Bamba et alj (12005b show that the 
enegy injected into accelerated electrons is always similar in 
young SN Rs, which is also c onsistent with our results. The- 
oretically, Berezhko & Volkl (12006) estimated the energy in- 
jected into accelerated protons to be 10^" ergs in the case of 
RX J1713-3946. This is also consistent with our results when 
we consider the electron and proton ratio of lO"'* in cosmic 
rays. 

At present, no direct information on the accelerated pro- 
tons is obtained by X-ray observations. One might consider 
whether we could estimate the total energy of accelerated pro- 
tons under the assumption of a similar A as for electrons. 
However, this is not possible because the injection efficiency 
is different in each case. On the other hand, the maximum en- 
ergy of accelerated protons, Ep"'"-\ can be estimated according 
to our model. The synchrotron and Tr^-decay loss timescales 
are negligible for protons, and thus Ep'""-^ is determined by the 
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condition face = fa^e, and we obtain 



(^aee) " 



4.8 X 10^ 



lO^c 



age 



10 Aioy V yr 



TeV 



(10) 

Although the effect of the wave damping and/or 
escape may be signifi cant i n older sy stems (e.g., 
IZirakashviU & Aharoniaiil 120071: lOhira et all 120101) . we 
ignore it here for simplicity. Figure |2] shows Ep"'"-^ as a 
function of fagg. In the free expansion phase, Ej,""" increases 
in proportion to the age, and protons are accelerated quickly 
up to around the knee energy in ^100 yr. When the shock 
velocity decreases, Ep'""" becomes smaller. Hence Ep'""-^ 
peaks at f^ge = f i . In contrast to the case of electron accel- 
eration, Ep"*"''{t\) (that is the maximum of Ep"""^) is larger 
for higher upstream density uq, and it changes only within a 
factor of ~ 2 for no = 0.1-5.0 cm"-'. Since oc no'''^ and 
t\ cx nQ~^l^, one finds Ep""'^(ti) oc B^ti oc no'''*- Protons are 
accelerated up to Ep""^(ti) ~ 10^^"'^ eV in any SNR, with 
little density dependence. This result might be important 
for explaining the break at ~ 10'^'^ eV in the observed 
cosmic -ray spectrum with the shock acceleration on SNRs. 



On the other hand, one can find E^' 



oc V 



,1/2 if 



Bel = Bamp V.5. Hence the maximum energy of electrons 
scales as Ee'""'' oc f° and oc f""*/'" for the free expansion and 
Sedov phases, respectively, so that it depends only weakly on 
the SNR age. 

Recent Fermi observations of middle-aged SNRs have 



shown the presence of a ga mma-ray spectral break around 
energies of a few GeV (e.g..lAbdoet"aLl l2009ll2010ah Sev- 
eral interpretations have been given (e.g ., 'Inoue et al. n2010t 
iLi & C hen 2010, 201 1; O hira. Murase. Y amazaki 201 1]), one 
of which is the escape of particles from acceleration sites. 
In this scenario, accelerated particles with energies of more 
than about 10-10 GeV have already escaped fro m the ac- 
celeration region (lOhira. Murase. Yamazaki 1 1201 lb . In con- 
trast, our present model has shown that even for fagg ^10"^ yr, 
~ 10 TeV electrons still exist around the shock front (see Fig- 
ure|2]l. This apparent discrepancy may come from the fact that 
most SNRs observed by Fermi are interacting with molecular 
clouds, whereas in our sample there is only one (W 28) under- 
going such a collision. Hence molecular clouds may play an 
important role in dispersing cosmic rays from the SNR shock 
into the interstellar matter. 
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Evolution of Synchrotron X-rays in SNRs 
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Fig . 1 . — Synchrotron X-ray luminosity in 2-10 keV band as a function of radius for Galactic SNRs. Solid lines represent the evolutionary model of synchrotron 
X-rays for no of 5.0 cm"^ (blue), 1.0 cm~' (red), and 0.1 cm~^ (green). We plot results only for the Sedov phase (<| < fagc < tj), in which the normalization of 
the electron energy distribution. A, is expected to be constant with time (see text for details). 
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Fig. 2. — The maximum energies of accelerated protons (solid) and electrons (dotted) as functions of the SNR age. The colors match those in Fig.[T] 
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TABLE 1 

Observational journals and physical parameters for SNRs with synchrotron 

X-RAYS. 



Target 


Distance 


Size 


Radius Lx" References 




(kpc) 


(arcmin) 


(pc) (lO^-* ergs s"') 



Gl. 9+0.3 

Cas A 

Tycho 

RXJ17 13-3946. 

SN 1006 

G28.6-0.1 

Vela Jr. 

W28 

CTB 37B 

G32.45+0.1 

G156.2+5.7 

RCW 86 



8.5 

3 4+0-3 

1.0 
2.2±0.1 

7 0+'-' 
0.65 ±0.35 
1.9±0.3 
10.2±3.5 
17^ 
1+0.3 
-0.7 

±0.4 



1.0 



1.5 
5 
8 

65 x55 

30 
13 X 9 

120 

48 

10"^ 
6 

110 
42 



1.9 

2.5*°-^ 

4 4+?:^ 
^•^-L3 
o 7+0.8 
^•'-0.7 

9.5 ±0.3 
11±2 
11±6 

13±2 
15±5 
15+5 



16!,, 
17±2 



4 g+O.l 
^•^-0.2 
110+f 
12+' 2 
6.5 

2 1+0.2 

0.88+|'| 
0.42^1 

0.0693;," 
97+3 ' 
4.5+"'-5 

060^-^^ ^ 

^•»-L0 



(1) 

(2) (3) (4) 
(5) (6) (7) 
(8) (9) 
(10) (11) (12) 
(13) 
(14) 
(15) (16) 
(17) (18) 
(19) 
(20) (21) (22) 
(23) (24) (25) 



NOT E. — (1) Reynolds et al. (2008); (2) Reed et al. : 1995); (3) Helder & Vink (2008); (4) Lazendic et al. 
(2001; (5)|Krause et al. ( 2008); (6) TamaRawa et al. ( 2009); (7) Seward et al. ( 1983); (8) Fukui et al. ( 2003); 
(9) Slane et al. ( 1999); (10) Winkler et al. : 2003); (II) Bamba et al. : 2008); (12) Yamasuchi et al. (200i^; (13) 
[Bambaet al. (2001); (14) Aharonian et al. (2007); (15) Velazquez et al. (2002); (16) Nakamura et al., sub mitted; 
(17) Caswell etal. (1975); (18) Nakamura et al. (2009); a9) |Yamag uchi et al. (2004); (20) Gerardv & Fesenl 
(2007); (21) Yamauchi et al. ( 1999); (22) Katsuda et al] 120091) ; (23) Rosado et al. : 1996); (24) Bamba et alj 
(200a) ; (25) i Vink etaL (2006) 
" In the 2-10 keV band. 

'' The upper-bound of X-ray luminosity is calculated with the assumption that the SNR uniformly emits syn- 
chrotron X-rays, since present observations do not cover the entire remnant. 
The size of radio partial shell is used. 



